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An updated version of an ITS-based phylogeny of the Calopteryginae, using 
sequences of 31 ingroup taxa, is given. The subfamily consists of 3 main clades, each 
with 2 subclades. Only clade 1 ( Calopteryx s. s.) is not exclusively Asian but extends 
to Europe and North America. In the East-Asian clade 2, the genus Matrona is found 
to be descended from an AtrocalopteryxAike ancestor. Several so-called South-East 
Asian Calopteryx probably either belong to Atrocalopteryx or to as yet unnamed 
genera near Atrocalopteryx. Archineura consists of 2 spp., limited to China and 
Indo-China, and is rather basal to clade 3. The subclade Neurobasis-Matronoides 
is worthy of further analysis. 


INTRODUCTION 

The dragonfly cohort Caloptera, composed of seven families, has recently been 
subjected to a phylogenetic analysis using the nuclear ribosomal genes 18S, 5.8S, 
and the internal transcribed spacers ITS1 and IRS2 (DUMONT et al., 2005). 
From this analysis, it appeared that the 5.8 and 18S genes in this group evolve 
at such a slow rate that almost all phylogenetically resolving power lies with the 
ITSses. 

The Calopterygidae and the Haeterinidae constitute the two most speciose 
families of the cohort, and within the Calopterygidae, the largest subunit is the 
Calopteryginae, with currently about nine genera recognized. 

In this subfamily, DUMONT et al. (2005) created a new genus, Atrocalopteryx , 
based on the position of a single species examined, A. atrata , distributed from the 
coastal fringe of Eastern Siberia to north-eastern China and Korea. However, this 
does not mean that the genus is monotypic: several species orginally described 
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in Calopteryx and occurring further south in China and Vietnam-Laos, are sus¬ 
pected to belong here, but none was so far available for molecular analysis. 

A Chinese endemic, Archineura incarnata, was found to cluster with Matrona 
and Atrocalopteryx, a surprising result, since by habitus and morphology, this 
taxon looked closely related to yet another calopterygine, “ Leucopteryx” ha- 
eterinoides, from Laos and Vietnam. The quotes suggest the uncertainty that 
surrounded (and still surrounds) the generic name Leucopteryx, since classi¬ 
cal morphology led to the expectation that haeterinoides was a true Archineura. 
Shortly after the publication of DUMONT et al. (2005), the suspicion about the 
placement of Archineura incarnata was confirmed, when it was found that this 
taxon had been identified and sequenced using a female that in fact belonged to 

Table I 

Alphabetical list of taxa and specimens examined 



Geographical origin 

Collector 

Archineura hetaerinoides 

Laos 

M. Hamalainen 

Archineura incarnata 

Nankunshan, Guangdong, China 

Boping Han 

Atrocalopteryx atrata 

Dokigawe River, Japan 

K. Inoue 

Atrocalopteryx atrocyana 

Nanling, Guangdong, China 

K.D.P. Wilson 

Atrocalopteryx atrocyana 

Tian Men Gao, Guangdong, China 

H.J. Dumont 

Calopteryx ’' coomani 

Vietnam 

M. Hamalainen 

Calopteryx aequabilis 

Wisconsin, USA 

S.W. Dunkle 

Calopteryx amata 

New Brunswick, Canada 

S.W. Dunkle 

Calopteryx Cornelia 

Yuragawe River, Japan 

K. Inoue 

Calopteryx exul 

Ifrane, Atlas, Marocco 

H.J. Dumont 

Calopteryx haemorrhoidalis 

Ifrane, Atlas, Morocco 

H.J. Dumont 

Calopteryx maculata 

Oklahoma, USA 

S.W Dunkle 

Calopteryx virgo virgo 

Laon, France 

H.J. Dumont 

Calopteryx xanthostoma 

Argens, Chateauvert, France 

M. Papazian 

Echo modesta 

Kanchanaburi, Thailand 

M. Hamalainen 

Matrona basilaris 

Omei Shan, Sichuan, China 

Su Rong 

Matrona basilaris 

Nankunshan, S. Guangdong, China 

H.J. Dumont 

Matrona basilaris 

Tian Men Gao, C. Guangdong, China 

H.J. Dumont 

Matrona ", basilaris ’* 

Hainan Island, China (2005) 

K.D.P Wilson 

Matrona ", basilaris ' 

Hainan Island, China (2006) 

H.J. Dumont 

Matrona cyanoptera 

Wash, Taipei, Taiwan 

W.C. Yeh 

Matrona nigripecta 

Chiang Mai, Thailand 

M. Hamalainen 

Matronoides cyaneipennis 

Mt. Kinabalu, Borneo 

M. Hamalainen 

Mnais andersoni 

Chiang Mai, Thailand 

M. Hamalainen 

Mnais mneme 

Hainan, China 

K.D.P Wilson 

Neurobasis chinensis 

Kanchanaburi, Thailand 

M. Hamalainen 

Neurobasis chinensis 

Asan lake, N. India 

H.J. Dumont 

Neurobasis chinensis 

Hainan Island, China 

K.D.P. Wilson 

Psolodesmus dorothea 

South Taiwan 

W.C. Yeh 

Psolodesmus mandarinus 

North Taiwan 

W.C Yeh 

Psolodesmus mandarinus kuroiwae 

Mt. Omoto, Ishikagi, Japan 

K. Watanabe 
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Matrona or A trocalopteryx. Thus, new material was collected from Guangdong, 
South China, to rectify this error. In the same general area (South China) and 
in Vietnam, more fresh specimens belonging to Matrona suspected to belong to 
Atrocalopteryx were collected as well, in order to substantiate better these subdi¬ 
visions of the calopterygine clade. 

Here, we sequence and analyse the internal transcribed spacers and interven¬ 
ing 5.8 S rDNA of these additional (seven in all) insects, with an aim at obtain¬ 
ing an improved phylogeny of the subfamily. Beside the Archineura , three were 
suspected to belong to Atrocalopteryx (two to the taxon atrocyana , and one to 
coomani ), and three to Matrona. They were inserted in a tree of Calopteryginae 
composed of 31 ingroup taxa, plus three outgroup taxa of the zygopteran fami¬ 
lies Chlorocyphidae, Diphlebiidae, and Megapodagrionidae. The origin and col¬ 
lectors of the specimens analysed is given in Table I. 

MATERIAL AND METHODS 

DNA EXTRACTION, PCR AMPLIFICATION AND SEQUENCING - The origins of the 
samples used in this study are listed in Table I. All were males that had been fixed and preserved in 
70-80% ethanol immeditaly upon collection ion the field. In the laboratory, muscular tissue was iso¬ 
lated from the synthorax and genomic DNA was extracted, using the protocol of the Puregene™ 
DNA isolation kit type D-5000A (Gentra Systems Inc., BlOzym, Netherlands),or using a modified 
CTAB protocol (KOCHER et al., 1989). Briefly, muscle tissue was crushed using a beadbeater and 
subsequently incubated for a minimum of 3 h at 60°C in 500 pi CTAB buffer to which 6 pi protei¬ 
nase K (10 mg/ml) was added. Next, 250 pi 7.5 M ammonium acetate was added and the mixture 
was centrifuged at 14,000 rpm for 10 min. The supernatant was saved and DNA was precipitated 
by adding an equal volume of isopropanol. The sediment was washed with 70% ethanol and redis¬ 
solved in 25 pi water. Small aliquots (usually 1 pi) were used as template for PCR. The complete re¬ 
gion separating the SSU and LSU genes and comprising the ribosomal spacers ITS1 and ITS2, and 
the conserved 5.8S gene, was amplified using the polymerase chain reaction. The total length of the 
fragment sequenced amounted to 712 bp. Finding primer pairs that yielded sequenceable amplicons 
was a tedious task. In all, we used the following primers: 

Vrain 2f : 5 -CTT TGT ACA CAC CGC CCG TCG CT- 3’ 

Ferris 2f : 5’-RGY AAA AGT CGT AAC AAG GT- 3’ 

Vrain 2r : 5’-TTT CAC TCG CCG TTA CTA AGG GAA TC- 3’ 

28R1 : 5’- TGA TAT GCT TAA NTT CAG CGG GT - 3’ 

5.8 n : 5’- TCG AAT TGT GAA CTG CAG GAC ACA T - 3’ 

5.8 r3 : 5’- TCC GTG GGC TGC AAT GTG CGT TCG AA - 3’ 

The PCR conditions were 30 s at 94°C, 30 s at 54°C and 2 min at 72°C for 40 cycles. After treat¬ 
ment with shrimp alkaline phosphatase and exonuclease I the PCR products could be used as tem¬ 
plate for cycle sequencing without further purification. 

PHYLOGENETIC ANALYSIS - The DNA sequences were aligned with MUSCLE (EDGAR, 
2004) using default settings. Analyses were performed under unweighted parsimony, maximum likeli¬ 
hood and Bayesian inference. Trees were displayed with TREEVIEW1.6.6 (PAGE, 1996). Parsimony 
analysis was performed using PAUP 4.0bl0 (SWOFFORD, 2003) with the following heuristic search 
settings: 100 random taxon addition replicates, followed by tree-bisection-reconnection branch (TBR) 
swapping (10 7 rearrangements). Gaps were treated as missing data. Nodal support was assessed by 
calculating bootstrap values from 100 bootstrap replicates obtained by heuristic search with 10 ran- 
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dom sequence additions each (FELSENSTEIN, 1985). 

Maximum likelihood analysis was performed using the general time-reversible substitution model 
with a gamma correction for among-site variation and corrected for invariable sites (GTR +1 + G). 
This model was identified as the best-fit model of DNA evolution for maximum likelihood analysis of 
the data-set by the likelihood-ratio test (LRT) and Akaike Information Criteria (AIC), implemented 
in ModelTest 3.7 (POSSADA & CRANDALL, 1998). Maximum likelihood was performed using 
the heuristic search option with stepwise taxon addition, TBR branch swapping, MulTrees option in 
effect, no steepest descent, and rearrangements limited to 100.000. Bootstrap values were determined 
from 100 bootstrap replicates by heuristic search with 10 random sequence addition replicates each 
(FELSENSTEIN, 1985). 

Bayesian analysis was performed using MrBayes, version 3.1.2 (HUELSENBECK & RONQUIST, 
2001). MrModeltest 2.2 (NYLANDER, 2004) also identified GTR + I + G as the best-fit model for 
Bayesian inference. The parameters for base frequencies, substitution rate matrix, gamma rate dis¬ 
tribution and shape and proportion of invariant sites were allowed to vary throughout the analysis. 
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oo 
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oo 
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Matrona cyanoptera (Taiwan) 

Matrona nignpecta (N Thailand) 
“Calopteryx ' coomani (North Vietnam) 
Atrocalopteryx atrata (Korea) 


r- Atrocalopteryx atrocyana (N Guangdong, China) 

{ 070 

Atrocalopteryx atrocyana (N Guangdong. China)] 
r- Neurobasis chinensis (Thailand) 

rl 040 

J ^Neurobasis chinensis (S. China) 


Neurobasis chinensis (India) 
Matronoides cyaneipennis (Borneo) 



Psolodesmus dorothea (Taiwan) 

097 

, Psolodesmus mandarinus kuroiwae (Japan) 

Psolodesmus mandarinus (Japan) 

Echo modesta (Thailand) 

Mnais andersoni (Thailand) 
oo 

Mnais mneme (Hainan, China) 

Archineura hetaerinoides (Laos) 
oo 

Archineura incamata (Guangdong, China) 
Philoganga vetusta (Hong Kong) 

56 

Heteragrion mitratum (Venezuela) 
Chlorocypha curia (Camerun) 


Cluster 1b 
N. American 
group 

Cluster la 
Temperate 
Eurasian 
group 


Cluster 2a 


Cluster 2 


Outgroup 


Fig. 1. Bayesian probability estimate of the phylogeny of the Calopteryginae based on internal tran¬ 
scribed spacer (ITS 1 and 2) data. 
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The Markov chain Monte Carlo process was run over four parallel chains (one cold and three heated) 
for 1,000,000 generations with trees being sampled every 100 generations. The bum-in value was set 
to 1000 trees (i.e. 100,000 generations) equating the number of generations needed to reach a stable 
value of all variable parameters in a preliminary run. Majority rule consensus trees were reconstruct¬ 
ed after discarding the bum-in. 


RESULTS 


We obtained the unambiguous sequences of the ITSes of all taxa listed in Ta¬ 
ble I, from which we computed four phylogenetic trees. The length of ITS varied 
froml63 to 203 bp, that of ITS2 from 206 to 229 bp, and that of 5.8 S from 163 to 
166 bp. The Bayesian and maximum-likelihood-based trees (heuristic tree only) are 
shown in Figs 1 and 2. Four most parsimonious trees were found, differing only 


rS 


Calopteryx aequabilis 
1000 

- 0 Calopteryx amata 
Calopteryx maculata 
Calopteryx exul 
Calopteryx xanthostoma 
Calopteryx haemorrhoidalis 
Calopteryx virgo virgo 


Calopteryx comelia 
Malrona basilaris (S. Guangdong) 
Matrona basilaris (N. Guangdong) 
Matrons basilaris (Sichuan) 

| |_ Matrona “basilaris” (Hainan) 

r^ u 1— Matrona “basilaris” (Hainan) 
Matrona cyanoptera 
Matrona nigripecta 


-e 


Atrocalopteryx atrocyana (Guangdong) 
1.5 

Atrocalopteryx atrocyana (Guangdong) 
Atrocalopteryx atrata 
“ Calopteryx “ coomani 



I— Neurobasis chinensis (Thailand) 
°Neurobasis chinensis (S. China) 

^- Neurobasis chinensis (India) 
Matronoides cyaneipennis 
Psolodesmus dorothea 

920 

psolodesmus mandarinus kuroiwae 
Psolodesmus mandarinus 
Echo modesta 
Mnais andersoni 

980 

Mnais mneme 

Archineura hetaerinoides 
860 

Archineura incarnate 

■ Chlorocypha curta 
- Heteragrion mitratum 

■ Philoganga vetusta 


Fig. 2. Maximum likelihood estimate of the phylogeny of the Calopteryginae based on internal 
transcribed spacer (ITS1 and 2) data. Bootstrap support based on 100 replicates and expressed as 
a percentage. 
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Fig. 3. Heuristic Maximum Parsimony tree with branch lengths of the 
phylogeny of the Calopteryginae based on internal transcribed spacer 
(ITS1 and 2) data. The MP analysis generated 4 most parsimonious trees 
of 1118 steps each. — [The scale bar represents 10 steps] 

Cluster 1 also has a subcluster that occurs in Canada and 


in minor details (the 
branching order of 
species within the 
genus Calopteryx ). 
Only the consensus 
tree is shown (Fig 
3). We also derived 
a neighbor-joining 
tree, but since this 
fully confirmed the 
previous ones, it 
is not reproduced 
here. 

The Caloptery¬ 
ginae form three 
well-defined clus¬ 
ters, each of which 
can be subdivided 
into two subclades. 
The vast majority 
occur in East and 
South-East Asia; 
only cluster 1 {Ca¬ 
lopteryx s. s.) ex¬ 
tends west across 
the southern half of 
Siberia and central 
Asia to reach Eu¬ 
rope and the Med¬ 
iterranean basin, 
the USA. 


DISCUSSION 


Most branches of our trees are well supported, allowing some rather robust 
conclusions to be drawn. Weaker support is only found in the relationship be¬ 
tween two Asian genera, Matrona and Atrocalopteryx, and some uncertainty also 
remains as to the descendence order within the genus Atrocalopteryx. 

All tree topologies reveal an undisputed monophyletic origin for the Caloptery¬ 
ginae (Bayesian probability 1.00; Maximum Parsimony, Maximum likelihood, 
and Neighbor joining all supported by 100% of the bootstraps executed), as es¬ 
tablished earlier by DUMONT et al. (2005), and the three main clusters found in 
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that paper are also confirmed (Figs 1-3). However, some important internal shifts 
fall to be recorded. Archineura incarnata is now closely clustering with " Leucop - 
teryx” haeterinoides , forming a subclade that is well supported (Bayesian prob¬ 
ability 0.99), as well as showing a close affinity between its two members (BP 
1.00). The generic status of Archineura is thus accepted, and the need for a sepa¬ 
rate genus 

(“Leucopteryx”) for taxon hetaerinoides lapses. New findings also include the fact 
that the position of Archineura is now fully embedded withinin cluster 3, and shows 
a sister-group relationship with the three other genera that consitute that cluster. 
Its relationship with such genera as Calopteryx, Matrona , and Atrocalopteryx is 
comparatively remote, not close as claimed in DUMONT et al. (2005). 

Cluster 2 is now better structured, with the clade Neurobasis-Matronoides in sis¬ 
ter position to a clade composed of Matrona and Atrocalopteryx. Matrona comes 
out of the analysis as a coherent, monophyletic genus, that is clearly descended 
from an Atrocalopteryx- like ancestor. The relationships within Atrocalopteryx are 
less straightforward, although clearly, “ Calopteryx ” coomani is seen to belong in 
Atrocalopteryx or a closely related taxon, and could well be the closest ancestor 
to Matrona. It certainly is not a true Calopteryx , and this position is supported 
by all three analyses independently. 

Cluster 2 is also characterized by only South-East Asiatic taxa, and suggests lo¬ 
cal evolution in the absence of major environmental disturbances such as glacia¬ 
tions. This involves quite a bit of cryptic speciation, well illustrated by Matrona , 
where the basilaris-g roup is clearly composed of several species that have so far 
not been well characterized morphologically (e.g. a taxon, provisionally called 
“basilaris ”, on the island of Hainan). In cluster 2b, the endemic Matronoides of 
Borneo is behaving as a genus distinct from Neurobasis , with the Indian speci¬ 
men of Neurobasis chinensis suggesting at least evolution at the subspecies level 
within this widespread oriental taxon. It will be interesting to add to the analysis 
representatives of Neurobasis from Borneo, Sulawesi and The Phillipines, where 
this genus is rather speciose, to further refine the relationships within this clus¬ 
ter. 

Cluster 2 furthermore raises the problem whether true Calopteryx , a genus of 
temperate and continental Eurasia and North America, really occurs south of the 
latitude of Japan. The little-known taxa C. melli, C. oberthueri and C. laosica from 
South China, Laos and Vietnam may well, like coomani , turn out to belong to 
different, still undescribed, genera. 

Cluster 1 ( Calopteryx ) and its two subclusters (an Eurasian and a North 
American one) remain unaltered as compared to the analysis of DUMONT et 
al. (2005), as well as the genera composing cluster 3 (except for the addition of 
Archineura). 
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CONCLUSION 

The subfamily Calopteryginae is composed of predominantly Asian repre¬ 
sentatives. A basal cluster of four genera (Psolodesmus, Echo, Mnais, and Arch- 
ineura) is currently restricted to subtropical-tropical East Asia. A second cluster 
(four genera, but probably one or two remaining to be defined) extends from the 
Asian tropics to temperate East Asia. Calopteryx , finally, is a monogeneric clus¬ 
ter that is relatively species-poor and the only component of the Calopteryginae 
that extends through North-Central Asia to Europe, reaching North Africa in 
the west, invading North America in the East. 
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